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Abstract 

We consider continuous time interlacements on Z d , d > 3, and investigate the 
scaling limit of their occupation times. In a suitable regime, referred to as the con- 
stant intensity regime, this brings Brownian interlacements on M. d into play, whereas 
in the high intensity regime the Gaussian free field shows up instead. We also inves- 
tigate the scaling limit of the isomorphism theorem of [40]. As a by-product, when 
d = 3, we obtain an isomorphism theorem for Brownian interlacements. 
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Introduction 



Informally, random interlacements provide a model for the local structure left at appro- 
priately chosen time scales by random walks on large recurrent graphs, which are locally 
transient. Their vacant set has non-trivial percolative properties, see for instance [9], 
[31], [38], [43], which are useful in the study of certain disconnection and fragmentation 
problems, see [6], [7], [37], [44]. Their connectivity properties have been actively investi- 
gated, see [5], [20], [28], [29], [30]. Random interlacements have further been helpful in 
some questions of cover times, see [1], [2]. They have also been linked to Poisson gases 
of Markovian loops, see [23], [41]. Recently, connections between random interlacements 
and Gaussian free fields have emerged, which underline the important role of occupation 
times of random interlacements, see [40]. 

In this article, we investigate the scaling limit of the field of occupation times of 
continuous time interlacements on Z d , d > 3. In the constant intensity regime, this brings 
into play the Brownian interlacements on M. d , d > 3. In the high intensity regime, the 
massless Gaussian free field shows up instead. Further, we investigate the scaling limit of 
the isomorphism theorem derived in [40] . Dimension three plays a special role, and when 
d = 3, we obtain as a limit an identity in law relating the occupation-time measure of 
Brownian interlacements on M 3 to the massless Gaussian free field, somewhat in the spirit 
of [22], [23], in the context of Poisson gases of Brownian loops at half- integer levels. 

We now discuss our results in more detail. We consider continuous time random 
interlacements on Z d , d > 3. In essence, this is a Poisson point process on a certain state 
space consisting of doubly infinite Z d -valued trajectories marked by their duration at each 
step, modulo time-shift. A non-negative parameter u plays the role of a multiplicative 
factor of the intensity measure of this Poisson point process, which is defined on a suitable 
canonical space, see [39], denoted here by (Q,A,P). The field of occupation times of 
random interlacements at level u is denoted by L x>u {uj), for x e Z d , m > 0, w 6 SI. It 
records the total duration spent at x by the trajectories modulo time-shift with label at 
most u in the cloud u, see [39]. 

We investigate the scaling limit of the field of occupation times of random interlace- 
ments on Z d , d > 3, and introduce the random measure on R d 

(0-1) CN = ^ £ L X>UN 5 § ,N>1, 

xez d 

where (un)n>i is a suitably chosen sequence of positive levels. In the constant intensity 
regime, that is when 

(0.2) u N = daN 2 ~ d , with a > 

(in this case the intensity measure Ef/^] of C N converges vaguely to ady as iV goes to 
infinity), we show in Theorem 3.2 that 

(0.3) C N converges in distribution to C a , as iV — > oo, 

where C a stands for the occupation-time measure of Brownian interlacements at level a, 
see (2.37). We construct Brownian interlacements in Section 2 using a similar strategy as 
in [38] or [42]. This is technically somewhat more involved in the present context, and the 
key identity is encapsulated in Lemma 2.1. The random measure C a is supported by the 
random closed set Xq , the Brownian fabric at level a, which is the union of the traces in 
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M. d of the doubly infinite trajectories modulo time-shift with label at most a in the Poisson 
cloud defining the Brownian interlacements, see (2.30). This random closed subset of M. d 
is a.s. connected when d — 3, but a.s. disconnected when d > 4, see Proposition 2.5. 

On the other hand in the high intensity regime, that is when 
(0.4) N d ~ 2 u N ->■ oo, 

we show in Theorem 3.3 and Corollary 3.5 that 



(0.5) C N = J 2N2 d du {£ N ~ E[£ N }) converges in distribution to the 

" UN massless Gaussian free field $ on M. d , 

which is the canonical generalized random field on S'(R d ), the space of tempered distri- 
butions on IR d , such that for any V in the Schwartz space iS(lR d ), the pairing ($, V) is a 
centered Gaussian variable with variance J V (y)G(y — y')V (y')dydy' , with G(-) the Green 
function of Brownian motion on IR d , d > 3, see (1.3). 

The low intensity regime, when N d ~ 2 UN — > 0, leads to a null limit for C N in (0.1), and 
will not be further discussed in the present work. 

From the isomorphism theorem for random interlacements, see [40], one knows that 
when ((p x )xez d is a discrete massless Gaussian free field, i.e. a centered Gaussian field with 
covariance E[(p x (p x r] = g(x — x'), with g(-) the Green function of simple random walk on 
Z d , see (1.1), independent from (L XtU ) xeZ d, one has for any u > the distributional 
identity: 

We explain in Section 4 how one quickly recovers (0.5) from this identity in the regime 
where, in essence, as N — > oo, the variance of Y1\ x \<n ¥x * s negligible compared to that 
of y/u N X]|x|<Af W This corresponds to the full range (0.4) when d — 3, but only the 
partial range of (0.4) where unN 2 log N — > oo, when d = 4, and unN 2 — >■ oo, when d> 5, 
see Remark 4.2. This fact singles out the special role of dimension 3, and motivates the 
investigation when d = 3 of the scaling limit of (0.6) (with adequate counter terms) in 
the constant intensity regime (0.2). Indeed, when d — 3, Theorem 5.1 roughly states that 
for any a > (denoting Wick products by : : , see Section 5) 

S -(v^ + v^%) 2; converges in distribution to :($ + v^a) 2 : , 

with mjv as in (0.2), and the last term defined by regularization of the massless Gaussian 
free field $ on M 3 (crucially using the fact that the Green function is locally square 
integrable when d — 3), see the beginning of Section 5. As a scaling limit of (0.6) we 
obtain the identity in law on 5'(R 3 ), for a > 0: 

(0.8) i:$ 2 :+£ Q = i:($ + v^) 2 :, 

where C a is independent of $ and denotes, as above, the occupation-time measure of 
Brownian interlacements on M 3 at level a, see Corollary 5.3. This identity in law has a 
similar spirit to some of the results in [22] , [23] in the context of Poisson gases of Brownian 
loops at half-integer levels, see Remark 5.4. 
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Let us describe how the article is organized. In Section 1 we collect some useful 
facts concerning Green functions, occupation times of random interlacements on Z d , and 
discrete as well as continuous free fields. Section 2 is of independent interest. It constructs 
the Brownian interlacements on M. d , d > 3, and derives some of their key properties, see 
Propositions 2.4, 2.5, 2.6. In Section 3, we investigate the limit of C N , respectively C N , 
in the constant intensity, respectively, high intensity regime. The main results appear 
in Theorems 3.2, 3.3 and Corollary 3.5. In the short Section 4 we use the isomorphism 
theorem (0.6) as a mean to recover (0.5). Section 5 focuses on the three-dimensional 
situation. The scaling limit of (0.6) (with adequate counter terms) is investigated in the 
constant intensity regime (0.2). The main results appear in Theorem 5.1 and Corollary 
5.3. 

Finally, let us explain our convention concerning constants. We denote by c,c',c,c 
positive constants changing from place to place, which simply depend on d. Numbered 
constants c , Ci, . . . refer to the value corresponding to their first appearance in the text. 
Dependence of constants on additional parameters appears in the notation. 

1 Some useful facts 

In this section, we introduce additional notation, and collect some useful facts concerning 
Green functions, occupation times of random interlacements on Z d , d > 3, and massless 
Gaussian free fields on Z d and R d , d > 3. 

We write | • |, respectively | • |oo, for the Euclidean, respectively, the supremum norm 
on R d . Throughout, we tacitly assume d > 3. We let B(y,r), y e R d , r > 0, stand for 
the closed Euclidean ball with center y and radius r. 

We denote by g(-, •) the Green function of simple random walk on Z d , that is, for 
x, x' G Z d , g(x,x') is the expected time spent at x' for the discrete time, simple random 
walk starting at x. The function g(-, •) is symmetric, and due to translation invariance, 
one has for x, x' in Z d 

(1.1) g(x,x') = g(x - x') = g(x' - x), where g(-) = g(-,0). 

One knows that g(-) < g(0) and that when x tends to infinity, cf. [21], p. 31, 

(1.2) g(x) ~ dG(x), where 

(1-3) G(y) = ^-L r(^ - l) \y\ 2 - d , for y e R d . 

Writing G(y,y') = G(y — y') = G(y' — y), one knows, see [10], p. 32, that the (0, oo]-valued 
function G(-, •) is the Green function of Brownian motion on R. d . 

For TV > 1, we denote by h N the lattice in R d : 
(1.4) h N = 1 Z d . 

For functions /, h on h N such that ^ ygLjv \f(y)h(y)\ < oo, we write 

<M> L =4 E f(y)Ky). 
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We rescale the Green function <?(•,•) and define 
(1.6) 9nM) = \ N d - 2 g(Ny,Ny'), for y,y' E L 



N, 



as well as g N {-) = g N (-, 0), so that g N (y, y') = g N {y - y') = g N {y' -y). By (1.2), (1.3), 
we also see that 

(1.7) lim sup \g N (y) - G(y)\ = 0, for every 7 > 0. 

N \y\>i 

We introduce the linear operator 

(1.8) G N f( y ) = * £ 9n(v, y') f(y'), y e W, 

Iy y'&L N 

which is well-defined when the function /: h N — y R is such that J^l^ ^(y, 1/') |/(y')l < 
00 for some (and hence all) y in Lat, in particular, when / vanishes outside a finite set. 

Similarly, when /, h are measurable functions on ~R d such that \fh\ is integrable, we 
write 

(1-9) (f,h)= [ f(y)h(y)dy. 

We also introduce the linear operator 

(1.10) Gf(y)= [ G(y,y')f(y')dy', for ye 



id 



which is well-defined when the measurable function / on R d satisfies J G(y, y') \f(y') \ dy' < 
00 for all y G M. d , in particular when / is bounded measurable and vanishes outside a 
bounded set. 

We now recall an identity for the Laplace transform of the field of occupation times 
{L x ,u) x & d °f continuous time random interlacements on 7h d at level u > 0. When 
V : Ljy — y K has finite support, the operator GjyV, which is the composition of the multi- 
plication by V with the operator Gjy in (1.8), sends bounded functions on L N into bounded 
functions on Ljv- We write HGw^IIl 00 -^ 00 f° r the corresponding operator norm (the space 
of bounded functions on Ljv being endowed with the sup-norm). One knows from Theo- 
rem 2.1 of [39] that when V: Ljv — > K has finite support and \\Gn \V\ \\l°°^>l°° < 1, then 
for u > 0, 

(1.11) l[exp{ E V(§) ^ 2 L X , U }] =e W {^N d -*(V,(I-G N V)-H) LN }. 



Note that when V vanishes outside a single point one readily finds by differentiation that: 
(1.12) I[£*,u] = u, for x e 7L d and u > 0. 

We then turn to the discussion of Gaussian free fields. Recall that we tacitly assume 
d > 3. We begin with the discrete case. We endow R zd with the product a-algebra, 
and denote by {(p x )xez d t ne canonical coordinates. The canonical law P 9 of the massless 
Gaussian free field on Z d is characterized by the fact that 



under P g , (ip x ) x & d is a centered Gaussian field with covariance 
E p9 [<p x <p x >\ = g(x,x'), for x,x' e Z d . 
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In the continuous case, we consider the space iS'(R d ) of tempered distributions on M. d . 
We endow ^'(R ') with the cylindrical a-algebra generated by the canonical pairings with 
functions of 5(R d ). We write $ for the canonical generalized random field on S'(M. d ) (i.e. 
the identity map). The symmetric bilinear form on <S(R d ) x l S(M a! ) defined by 

(1.14) E(V,W) = J V(y)G(y-y')W(y')dydy'{ = (V,GW) = (GV,W)) 

is positive definite (see for instance [36], p. 75). It satisfies the bound 

\E(V, W)\<c \\V\\ LHRd) (\\W\\ Lx{Rd) + H^IUi^)), for V,We S(R d ), 

and hence is continuous (we endow iS(M a! ) with its usual Frechet topology, see [17], p. 6-8). 
By Minlos' Theorem, see Theorem 2.3, p. 12 of [34], or Theorem 2.4.1, p. 28 of [17], there 
exists a unique probability measure P G on S'(R d ) such that 

under P G , for each V G S(R d ), <$, V) is a centered Gaussian variable 
with variance E(V, V). 

The law P G describes the massless Gaussian free field on R d , see Chapter 6 §2 of [15]. 



2 Brownian interlacements 

In this section we construct the Brownian interlacements on R d , d > 3, introduce their 
occupation-time measure, and discuss some key properties of these objects. We follow 
a similar approach as in [38], [42]. However, the continuous set-up is technically more 
challenging. An important step is the construction of the intensity measure of our basic 
point process on a space of doubly infinite trajectories modulo time-shift. We mainly rely 
on a crucial compatibility property of the measure expressed in "local charts" , see Lemma 
2.1 (cf. Theorem 1.1 of [38] and Theorem 2.1 of [42]), rather than on an approach based 
on projective limits, see [45], [32], following the outline of [16]. 

We first need to introduce notation. We recall that we tacitly assume d > 3. We denote 
by W+ the subspace of C(R + ,R d ) of continuous Revalued trajectories tending to infinity 
at infinite times, and by W the subspace of C (R, R d ) consisting of continuous bilateral 
trajectories from R into M. d , which tend to infinity at plus and minus infinite times. We 
write X t , t > 0, and X t , t 6 E, for the respective canonical processes, and denote by 6 t , 
t > 0, and t , t G R, the respective canonical shifts. The spaces W + and W are endowed 
with the respective a-algebras W + and W generated by the canonical processes. Given 
an open subset U of lR d , w G W+, we write Tjjiw) = inf{s > 0;X s (w) ^ U} for the exit 
time of U. When F is a closed subset of M. d , we write Hp{w) = inf {s > 0; X s (w) G F)}, 
H F {w) = inf{s > 0;X s (w) G F}, for the respective entrance time and hitting time of 
F. When w G W, we define Hp{w) and Tu{w) by similar formulas, simply replacing the 
condition s > 0, by s G R. 

We then consider W* the set of equivalence classes of trajectories in W modulo time- 
shift, i.e. W* = W I ~, where w ~ w', when w(-) = w'(- + t) for some t G R. We denote by 
7r* the canonical map W — > W*, and introduce the a-algebra W* = {A C W*; (7r*) _1 (A) G 

W}, which is the largest a-algebra on W* such that (W, W) (W*, W*) is measurable. 
Incidentally, note that a random variable Z on W invariant under 9 (i.e. Z o 6 t = Z, for 
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all t in M) determines a unique random variable Z* on W* such that Z* o tt* = Z. Given 
a compact subset K of ~R d , we write Wk for the subset of trajectories of W that enter K, 
and Wk for its image under 7r*. 

Since d > 3, and Brownian motion on IR d is transient, we view P y , the Wiener mea- 
sure starting from y G M. d , as defined on (W+, W+), and write ^ for the corresponding 
expectation. When p is a finite measure on R d , we write P p for the Wiener measure with 
"initial distribution" p and E p for the corresponding expectation. 

When B is a closed Euclidean ball of positive radius (as a shorthand in what follows, we 
will simply write that B is a closed ball), and when y ^ B, we define -P<f [•] — • I -5b — 
oo] to be the law of Brownian motion starting at y conditioned never to enter B, and 
write i?^ for the corresponding expectation. When y G (95 (the boundary of B), as 2; 
in R d \B tends to y, the measures Pf converge weakly (on C(IR + ,IR d )) to a measure P® 
supported on {w G W + ;w(0) = y and w(t) ^ i? for all i > 0}, which can be represented 
as the Brownian excursion measure in M. d \B starting from y conditioned on the event of 
finite positive mass that {Hb = 00} (see Theorems 4.1 and 2.2 of [4], and also Theorem 
3.1 of [8] for bounds on the exit time from small balls centered at y under Pf). Using 
rotational invariance and the explicit nature of the conditioning, one can also in a more 
direct fashion establish the above mentioned facts using a skew product representation of 
the law Pf. 

We write pt{y,y') = ^ 2 Jt) d / 2 ex P{~ ^~2t ' }' ^ > 0, y,y' G M. d , for the Brownian 

transition density, so that with similar notation as below (1.3) 

POO 

(2.1) G(y,y')= / p t (y,y')dt, fory,y'eR d . 

Jo 

Given a compact subset K of R d , we denote by the equilibrium measure of K (see 
Proposition 3.3, p. 58 of [36], or Theorem 1.10, p. 58 of [27]). This measure does not 
charge points of K and is supported by the boundary of K (actually, the boundary of 
the unbounded component of K c ). Its total mass is the capacity cap(i^) of K. The 
equilibrium measure ex is related to the time of the last visit of K via (see Theorem 3.4, 
p. 60 of [36]): 

(2.2) P y [L K >0,L K e dt, X Lk G dz\ = p t (y, z) e K {dz) dt, 

where for w G W + we write L K {w) = sup{t > 0;w(t) G K} (with the convention 
sup = 0), for the time of last visit of w to K. 

For B a closed ball (see above (2.1) for the terminology), we introduce the following 
(finite) measure on W% = {w G W;H B {w) = 0} (C Wb), the subset of W of bilateral 
trajectories that enter B at time for the first time, 

(2.3) Q B [(X_ t ) t > G A', X G dy, (X t ) t > G A] = e B (dy) P y B [A'] P y [A], for A, A' G W+. 

The next lemma states an important compatibility property of the above collection of 
finite measures. In essence, the measures Q B should be thought of as the expressions in 
the (natural) "local charts" of W* with respective domains W B (C W B ) of the intensity 
measure we are trying to construct (see also below (1.15) of [38]). 

Lemma 2.1. Assume that B,B' are closed balls, and B lies in the interior of B' , then 

(2.4) 6 Hb o(l{H B < 00} Q B ,) = Q B . 
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We postpone the proof of Lemma 2.1 for the time being, and explain how the con- 
struction of the intensity measure of our basic Poisson point process proceeds. We first 
observe that for K compact subset of R d we can unambiguously define 

(2.5) Q K = 6 Hk o (1{H k < 00} Q B ), for any closed ball BDK, 

and this definition agrees with the definition in (2.3) when K is a closed ball. Indeed, 
when the closed balls B\ and B 2 contain K, we can find a closed ball B' containing BiUB 2 
in its interior, so that by Lemma 2.1 we have 

9 Hk o (1{H k < oo}Q Bl ) = 9 Hk (1{H k < oo}9 Hbi o (1{H Bi < 00} Qb 1 )) 

(2.6) = 6 Hk o {9 Hbi (1{H k < 00} Q B ,)) 

= 9 Hk o (1{H k < 00} Q B ,), 

and a similar identity holds with B 2 in place of B\. Hence (2.5) is unambiguous and 
agrees with (2.3) when K is a closed ball. 

We now come to the theorem constructing what will be in essence up to a non- 
negative multiplicative factor the intensity measure of the Poisson point process defining 
the Brownian interlacements. 

Theorem 2.2. There exists a unique a- finite measure v on (W*, W*) such that for each 
compact subset K of M, d 

(2.7) l Wk u = 7i*oQ K . 

Proof. Given a sequence of compact subsets K n f R d , we have W* = \J n>0 W Kn , and the 
uniqueness of v is immediate. For the existence of is, it suffices to check Fhat for K C K' 
compact subsets of M d 

(2.8) ti*oQ k = tt*o(1 Wk Q k ,). 

Indeed, one then defines v so that lw* v = n* o Q K for some K n f noting that the 
restriction to W Kn of rr* oQ Kn+1 equals ir* o (l\v Kn QK n+1 ) which equals tv* oQ Kn by (2.8). 
Using (2.8) once again, one sees that v does not depend on the sequence K n) and (2.7) 
immediately follows (choosing K = K). To prove (2.8), we note that when the closed 
ball B contains K', then by (2.5) we have 

Qk = 9h k (IwkQb) and 1w k Qk' = ^w k ®h k , (\w k ,Qb) = 0h k , (^w k Qb)- 
Hence, the images under n* of Qk and 1w' k Qk' coincide, that is (2.7) holds. □ 

We are now reduced to the 
Proof of Lemma 2. 1: 

It suffices to show that for any continuous compactly supported function /:]&—>■ M d , 
j e B/ (dy')E y ,[H B < 00, ^-H B yx vdv ^ ^j- /(„).*„*] ^ Hb) ] = 

(2-9) J eB (dy) E* [e 1 ^ /HO-*-*] E y [ e */o~ /<«)■*.*] , whe re 
F'(y', t) = E^' [e^o°° /(-«-t)-x««fo] ? for y , e dB ^ f > q 
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Indeed the second line of (2.9) coincides with E^ B [e l ^^ v ^' Xvdv ], whereas by the strong 
Markov property at time H B , the first line of (2.9) equals 

E Qb '[H b < oO j e i ^f( V ~ HB " > ' Xvdv ] = E 9H B 0< - 1 { HB<00 } Q Bi)[ e i fmf( v )- x vdv^ 

and the claim (2.4) follows. We will thus prove (2.9). We denote by P* , the Brownian 
bridge measure in time t > 0, from y to y' in M, d , see [36], p. 137-140, and write E yyl for 
the corresponding expectation. We consider y e B' and use the last exit decomposition 
at the last visit of B', see Theorem 2.12 of [14] and (3.39) on p. 60 of [36], to assert that 
under P y , conditionally on Lb 1 — t (with t > 0), and X Lgl = y' (with y' G dB), the law 
of (X Lb , +v ) v > and (X v ) < v < t are independent, respectively distributed as P yl and P yy , 
and (Lb>, X Lb ,) has distribution pt(y, y') eB'(dy') l{t > 0} dt. Writing U = B c we find as 
a result that for y in B' 

E y [e i fS°X- v > Xvdv ,HB = oo] = 

dtj e^)^'"'^"^ 
jTcft J eB,(dy')F'(y>,t)E t y , !y [eti^- t > x " d \T u >t] p t (y',y), 

where, in the last step, we used that the law of (JT t _.) <.<t under P* , equals P y , , 
cf. (A. 12), p. 139 of [36], that pt(-, •) is symmetric, and the notation from below (2.9). 

We introduce the notation F(z) = ip(z) E z [e 1 ^ fH-x v dv^ f or z e jg>^ w ]j ere ^ j s a con _ 
tinuous compactly supported [0, l]-valued function, which equals to 1 on a neighborhood 
of P. We pick e > and introduce the finite measure 

m £ (dy) = - dz F(z) P Z [H B <e,X £ e dy] 

(2.H) \ \ 

= t J dz F(z) P Z JH B < e] p £ {z, y) dy. 

We integrate (2.10) with respect to m £ (dy). The first line of (2.10) yields 

J dzF(z) \ E Z [H B < e, E Xe [H B = oo,e^/(-^]] 

/ dzF(z) ± E Z [H B < e, P Xs [H B = oo] E* [ e *r /HO-*-*]] 

J £ ' - J property 

J dzF(z) \ E Z [H B <e,H B o6 £ = oo, EB t [e i f"H- v > x " dv ]] = 
J dzF(z) - £ E z [0<L B < e, PfJe^o°° . 

By the last exit decomposition of Brownian motion starting at z at the last visit of P, 
the above expression equals 

dzF(z) 1 - f dsp s (z,y)e B (dy)E*[El ySS°W.*>]] = 
6 Jo 

e B {dy)- £ jT<fe( J p s (y,z)F(z)dz)E^[ElJe i ^f(-y x ^]] (by symmetry of p s (; ■)). 
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Using the fact that for y G dB, P^-a.s., as u — > 0, Pj^ converges weakly to P^ (one 
could actually also invoke the Markov property under P® at this point), and dominated 
convergence, as e — > 0, the above quantity tends to 

(2.12) J e B (dy) E y [e* /o°° ^^*] Pf [e^ 

Using the symmetry of p e {-, •) and the fact that for y G U = B c , P e zy \E B < s] — P y , z [H B < 

e] (when z G B, note that P z -a.s., H B = 0, and both terms equal 1), the last line of (2.10) 
after integration with respect to m £ (dy) yields 

J e B ,(dy>) J™^F'(y',t) J dzdyF(z) 

E t',v ¥ !t ° f{v - tyXvdv , Tu > t] Pt (y', y) P^[Hb < e] p £ (y, z) = 
J e BI {d^J^^F t {^t)E^[e i ^^ t > x ^ 1 t<H B <t + e 1 F(X t+e )} = 

[ e B ,(dy')E y ,\H B < oo, - ["* e 1 ^ ^ v ~ tyXvdv F(X t+£ ) F'(y', t) dt 

J 1 £ J(H R -e) + 



-H B 
'(H B -s) 

Applying dominated convergence as e — > 0, the above quantity tends to 



(2.13) J e B ,{dy')E yl [H B < oo, e^ B ^E^^ m*«**\ F '( y ', Hb )] . 

Comparing (2.12) and (2.13), we have shown (2.9). This proves Lemma 2.1. □ 
Remark 2.3. 

1) When B is a closed ball, for w G W B we can define L B (w) = sup{t G M.,X t (w) G B}. 
Applying the last exit decomposition formula, as in (2.10), to the forward term in (2.3), 
we see that for A,A'<E W + we have 

Q B [(X- t ) t >o G A', (X s )o< s < Lb G •, (X Lb+s ) s > G A] = 
(2 ' 14) dt J e B (dy') e B (dy) P*[A'\ P^ y [\ P y B [A] p t (y', y) 

(where (X s ) < s <l b is viewed as a random element of the space T of continuous Revalued 
trajectories of positive finite duration, measurably identified with C([0, 1], M. d ) x (0, oo) 
via the map (w,T) — > w(^) G T). 

From this identity using the symmetry of p t (-, ■) and the already used fact that P y , 
is the image of P y y , under time reversal, we find (analogously to (1.40) of [38]) that 

(2.15) (A Ls _() ieR under Q B has law Q B . 

One can introduce on W* the time inversion w* — > w*, where w* = n*(w), with w(t) = 
w(—t) for t 6 1, and w G W arbitrary such that tc*(w) = w*. It follows from (2.15) 
that the image (tt* o Q B ) y of (ir* o Q B ) under time inversion coincides with tc* o Q B . By 

Theorem 2.2 it follows that the image v of v under time inversion satisfies 

(2.16) v = v. 



9 



2) In an easier fashion one sees that when one considers y G M d and denotes by r y the 
translation by — y on W*, which to w* associates w* — y (defined in an obvious way), or 
a linear isometry R of M. d , and its natural operation w* — > Rw* on W*, we find that 

(2.17) v is invariant under r y , 

(2.18) v is invariant under R, 

3) When A > 0, one can define the scaling s\(w) = Aw(p-) for w G W, and on W* via 
S\(w*) = 7i*(sx(w)), for w E W arbitrary such that ir*(w) = w*. From the identity valid 
for arbitrary A > and closed ball B 

(2.19) s x oQ B = \ 2 - d Q XB , 
one finds by Theorem 2.2 that 

(2.20) s x o v = \ 2 ~ d u, for A > 0. 

4) One knows that when the compact set K(0 M. d ) lies in the interior of the closed ball 
B, then e K (-) = P eg [H K < oo,X Hk G •], see Theorem 1.10, p. 58 of [27]. By (2.3) and 
(2.5) we thus see that 

(2.21) Q K [(X t ) t > G •] = P eK . 

We can now introduce the measurable map from into W+ defined by w* G — > 
w *,k,+ = ( w (h k + t)) t>Q , for any w e W K with tt*(w) = w*. Then, by (2.7) and (2.21) 
we see that 

(2.22) the image of lw^v under w* — > w*' K ' + equals P eK . 

□ 

As a next step we introduce the canonical space for the Brownian interlacement point 
process, namely the space of point measures on W* x M + , 



(2.23) VI = {lo = Y, fyof.ai), wi th «, on) G W* x [0, oo) and u{W* K x [0, a}) < oo, 
for any compact subset K of R d and a > } . 



We endow f2 with the cr-algebra A generated by the evaluation maps w — > ui(B), for 
B G W* ® B(M. + ), and denote by P the law on (£l,A) of the Poisson point measure with 
intensity measure v®da on W* x M + . For K compact subset of M. d and a > 0, we consider 
the measurable function on Q with value in the set of finite point measures on W + 

(2.24) HK, a ( u ) = E H«i < «>< e W* K } 5 whK:+ , when w = ^ G 

i>0 i>0 

with the notation from above (2.22), We readily see by (2.22) that 

(2.25) jj,K,a is a Poisson point process on W + with intensity measure aP eK . 

As a direct consequence of (2.16) - (2.20) we obtain the following invariance properties. 
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Proposition 2.4. P is invariant under the following transformations on fi: 



(2 26) 


1 1 1 — 

wj — 


u (w*,ai) ' 

i>0 








OJ — 


- X. . * 
u (w*,a l ) 7 

i>0 


u («)*-2/,ai) 

j>0 


( WiblL y t JK. y, 


(2.28) 


OJ - 


~- E — ' 

i>0 


' E d(Rw*, ai ) 
i>0 


fwitt i? linear isometry ofW 1 ), 


(2.29) 


OJ - 


i>0 


■ E 5 ( S a«),A2- 
i>0 




When oj 


en, a 


> 0, r > 0, the formula (see 


the beginning of Section 1 for notation) 


(2.30) 




2?H = U 

i>0:«i 

and 7r*(wj) = 


U B(wi(s),r), where w G £ 

<a seM i>0 

w* for i > 0, 



defines the Brownian interlacement at level a with radius r. By construction, see (2.23), 
this is a closed subset of ~R d . When r = 0, we refer to Xq (cu) as the Brownian fabric at level 
a. The terminology is truly pertinent when d — 3, see (2.36) below. The complement of 
X"(w) is an open subset of M. d , the vacant set at level a and radius r: 



(2.31) V?(oj) = R d \l?(oj), for oj e n, a > 0, r > 0. 

The set E of closed (possibly empty) subsets of R d can be endowed with the cr-algebra T 
generated by the sets {F e E; F n If = 0}, where if varies over the compact subsets of 
]R d , see [25], p. 27 (incidentally, this is the Borel a-algebra for a metrizeable topology for 
which E is compact, see [25], p. 3 and 27). The E-valued map X" on f2 is measurable, 
and one can thus consider its law on (E, J 7 ). 

Proposition 2.5. (ce > ; r > ; y G M d ; i? Zmear isometry, \ > 0) 
Qr is determined by the identity 

(2.32) Q"({F G E; F n K = 0}) = P[X r a n K = 0] = e - Qcap ^, for K CR d compact, 

where K r = K + 5(0, r) £/ie se£ of points within \ • \ -distance at most r from K. Moreover, 
under P, 

(2.33) X" + y nas same /aw as X" (translation invariance), 

(2.34) R(%r) has same law as X" (isotropy), 

(2.35) AX" aas same law as X^ r da (scaling), 

(2.36) Xq is a. s.- connected, when d = 3, and disconnected, when d > 4 and a > 0. 
Proo/. Note that {X r a n K = 0} = {oj G x [0, a}) = 0}, so that 

P[X r Q n K = 0] = c- "^) (2 - 7) i?- 21) e -«ca P (jf r ) > whence (2.32). 
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The identity (2.32) determines on a 7r-system generating J 7 , and the first claim follows. 
The identities (2.33) - (2.35) are direct consequences of (2.27) - (2.29) (alternatively, they 
can be derived from (2.32)). The last claim is a direct consequence of the fact that 
for any closed ball B,Q B (dw) <g> QB(dw')-a.s., w(M) D w'(M) ^ 0, when d — 3, but 
w(M) fl = 0, when d > 4, as a consequence of the fact that Brownian paths meet 

each other when d — 3, and, when starting from different points, miss each other when 
d > 4, see [11] (when d > 4 we also use (2.15) to take care of the bilateral nature of the 
paths). □ 

As a last topic of this section, we introduce the occupation-time measure of Brownian 
interlacements at level a > and discuss some of its properties. It is the locally finite (or 
Radon) measure on R d defined for u e Q, a > 0, A e B(R d ) by 

C a (u)(A) = / H w i( s ) £ A}ds, where u = Y ^«,<*;) witn 

(O 17\ i>0:ai<a JR */ \ * r ■ ^ n i>0 

l^-J'J " 7T (lUj) = IU*, for I > 0, 

= (u,fA® l[0,a]>, 

where }a{w*) = f R lA(w(s))ds, for mj 6 If arbitrary such that 7r*(u>) = w*, and (u>,h) 
stands for the integral of h with respect to uo. Note that by the choice of Q and W*, the 
expression in (2.37) is finite when A is a bounded set. From the second line of (2.37), the 
dependence in u is ^4-measurable, and C a defines a random measure on IR d , see Chapter 
1 of [19]. Further when A e B(M. d ) is bounded and B a closed ball containing A, then for 
a > 

E[£ a (u>)] = E[(u>, f A ® l [0) a]>] = a(^, /a) 

(2.38) (2 - 7 U 2 - 3) aE eB [J™l A (X s )ds\ = aj e B (dy) G(y,y') l A (y') dy' 

— a \A\ (with \A\ the Lebesgue measure of A), 

since / e B (dy) G(y,y') = 1 for y' E B, by (2.2), (2.1). In other words: 

(2.39) the intensity measure of C a equals ady. 
Observe also that by inspection of (2.37) 

(2.40) the support of C a coincides with Xq . 

The next result will be very useful and provides an expression for the Laplace transform 
of C a . 

Proposition 2.6. When V is a bounded, measurable, compactly supported function on 
R d , and 

(2.41) || C | V"| || L °°(Rd) < 1 (see (1.10) for notation), 

then I — GV is a bounded invertible operator on L°°(M. d ) and for any a > 0, 

(2.42) E[exp{(C a ,V)}} = exp {a(V, (I - GVy 1 !)}, 

with the notation (1.9) and (C a ,V) the integral ofV with respect to C a . 
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Proof. Up to routine modifications, this follows by similar arguments as Theorem 2.1 of 
[39]. □ 

Remark 2.7. As a direct application of the invariance property (2.29) and the definition 
(2.9) we see that for a > 0, the occupation-time measure C a has the following scaling 
invariance 

(2.43) C a '= A 2 h\ o C xd -2 a , for A > 0, 

where h\(y) = \y denotes the homothety of ratio A on M. d . One can also recover the 
identity (2.43) from (2.42). □ 



3 Scaling limits of occupation times 

In this section we study the scaling limit of the random field of occupation times of 
random interlacements on Z d , d > 3. In the constant intensity regime (0.2), we show 
that C N (see (0.1)) converges in distribution to the occupation-time measure of Brownian 
interlacements, see Theorem 3.2. In the high intensity regime (0.4), we show that C N (see 
(0.5)) converges in distribution to the massless Gaussian free field, see Theorem 3.3 and 
Corollary 3.5. 

We tacitly endow the set of Radon measures on R d with the topology of the vague 
convergence (see A. 7 of [19]). Given a positive sequence (un)n>i, it follows from (1.12) 
that 

(3.1) E[C N ([0,l) d )] = 1 1 N d - 2 u N . 

This quantity equals a(> 0) in the constant intensity regime (0.2), and tends to oo in 
the high intensity regime (0.4). The low intensity regime (when N d ~ 2 UN —> 0) leads to a 
convergence in distribution of C N to the null measure on M. d , and we will only focus on 
the constant and high intensity regimes. 

We begin with some preparation. We consider an integer M > 1, and 

(3.2) a continuous function V on R d with support in C M = [-M, M] d . 

We recall the notation (1.8), (1.10) from Section 1. With a slight abuse of notation, we 
still denote by V the restriction of V to (see (1.4)). 

Proposition 3.1. (under (3.2)) 

(3.3) sup \\G N \V\ ||l~ (Ljv ) < c (M) \\V\\ L ~ { c M ni. N ) < c {M) \\V\\ L oo (Rd) . 

N>1 

(3.4) hm<\/,(G iV \/r- 1 l> Ljv = <\/,(G\/r- 1 l>, for all n > 1. 

Proof. We begin with the proof of (3.3), and assume, without loss of generality, that 

IIUII 

L°°(c M riL N ) — 1- By a classical harmonicity (or martingale) argument one sees that 
^Pj/eLjv \ G N\V\(y) \ = sup„ 6C - Mr(Ljv \G N \V\(y)\, since V is supported in C M . It then follows 
from (1.2), (1.8) that 

sup \G N \V\(y)\<±(l+ E ]^)< C ' M2 

(3.5) yec M nh N iV ^ 0<|x| oo <2MAr Fl ' 

(x belongs to Z d in the sum), 
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and the claim (3.3) follows. Note that the continuity of V was not needed for the proof 
of (3.3). We then turn to the proof of (3.4). As we now explain, it suffices to show that 
for all bounded continuous functions W on M. d , 

(3.6) lim sup \(G N VW)(y) - (GVW)(y)\ = 0. 

N yeC M riL N 

Indeed, once (3.6) is established, we note that GVW is bounded continuous on M. d as 
a convolution of the locally integrable function G(-) (see (1.3)) with the compactly sup- 
ported continuous function VW. Then, we observe that for n > 1, 

sup \(G N V) n l(y) - (GV) n l(y)\ < sup ((^^((G^VT" 1 ! - (CV^ 1 ^^ 

yeC M nh N yeC M nh N 

sup KGxViGVr- 1 !)^) - (GViGVr-H^y)]. 

yeC M nh N 

Using induction over n > 1, we see that the last term tends to zero with N, by (3.6) 
and the observation below (3.6), and the first term after the inequality sign tends to zero 
with N by the induction hypothesis and the first inequality in (3.3). Thus, once (3.6) is 
proved, it follows that 

(3.7) sup |(G iV \/) n - 1 l(-) - (G\/) n - 1 l(-)| — >■ 0, for each n > 1. 

C M nh N n 

The claim (3.4) now follows by a straightforward Riemann sum approximation. 

We now prove (3.6) and introduce the shorthand F = VW. The function F is con- 
tinuous on R d with support in C M - Then, for 7 G (0, 1), N > 1, and y G C M H h N , we 
have 



< 



\G N F(y)-GF(y)\= -L £ g N ( y >) F(y - y') - [ G(y) F(y - y') dy' 

JV y'€L N JR d 

(3.8) E 9n{v')+ [ G(y')dy')\\F\\ L ^ (Rd) + 

4 E 9N(y') F(y -y')- f G(y')F(y - y') dy' . 

7<l?/'|oo<2M ^7<|J/'| 00 <2M 

The functions G(-) F(y — ■) are uniformly continuous on {y' G M. d ; 7 < \y'\oo < 2M} 
uniformly in y G Cm- By (1.7) and a Riemann sum approximation argument we see that 
the last term of (3.8) tends to zero uniformly in y G Cm H Ljy. On the other hand, the 
second line of (3.8) is bounded by (-^(^N 2 + 1) + c'7 2 ) ||-F||.L°c,( M d). Taking a limsup in 
iV, and then letting 7 tend to zero, we obtain (3.6). The claim (3.4) now follows, and this 
completes the proof of Proposition 3.1. □ 

We are now ready to state the convergence result for the law of C N in the constant 
intensity regime (i.e. N d ~ 2 UM = da, see (0.2)). The definition of C a appears in (2.37). 

Theorem 3.2. (under (0.2)) 

(3.9) C N converges in distribution to C a , as N — > 00. 
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Proof. By Theorem 4.2, p. 22 of [19], it suffices to show that for any continuous compactly 
supported function V on M, d 

(3.10) (C N , V) converges in distribution to (C a , V), as iV — > oo. 

Without loss of generality we assume that V satisfies (3.2) and ||V|| LO o( R d) < 1. By (1.11) 
and (3.3) we see that for \z\ < ^ (with c from (3.3)) and N > 1, 

E[exp{z(C N , V)}} = exp {a(zV, (/ - zG N V)- l l)^} 
(3 ' U) =e W {aJ2z-(V,(G N Vr-'l) hN }. 

n>l 

In particular this shows (with (3.3)) that 

(3.12) sup E[ cosh (r(C N ,V))} < oo, when r < Cq \ 

N>1 

Thus, the laws of (C N ,V), N > 1, are tight, and the random variables e z ^ N ' v \ N > 1, 
\Rez\ < r(< Cq 1 ), are uniformly integrable. Hence, if along some subsequence Nk, k > 1, 
the random variables (C N , V) converge in distribution to a random variable U , it follows 
from Theorem 5.4, p. 32 of [3], that for \z\ < c 1 in C one has 

E[e zU ] = lim E[e*< £jVfe ' y >] 

k 

(3.13) (3 = 1} >r ex p { a ^ zn ^ ^^) n " 1 > Ljv } 

(3 - 3) i 3 - 4) exp {a E z n (V, {GVy- 1 !) } ( = 2) E[e z ^}, 

n>l J 

where in the last step we have used the fact that ||G \V\ \\l°°(k<i) < c as a result of (3.3) 
and (3.6). By analyticity, the first and the last expression in (3.13) are equal in the strip 
\Rez\ < c^ 1 of the complex plane. Hence U and (C a ,V) have the same characteristic 
function, and we have shown (3.10). This completes the proof of Theorem 3.2. □ 

We now turn to the discussion of the high intensity regime (0.4) (i.e. N d ~ 2 u^ — > oo), 
and recall the definition of the random signed measure C N from (0.5): 



For a bounded measurable function V on M. d with compact support, the notation (C N , V) 
will refer to the integral of V with respect to C N . In essence, the next result states that 
C N tends in distribution to the massless Gaussian free field, see (1.15). This fact will be 
made precise in Corollary 3.5 below. 

Theorem 3.3. (under (0.4)) 

When V is continuous compactly supported function on M. d , then, as N — )• oo ; 

(3.15) (C N ,V) converges in distribution to a centered Gaussian variable 
with variance E{V, V) (1 = } J V(y) G(y - y') V(y') dydy'. 
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Proof. Without loss of generality, we assume, for convenience, that (3.2) holds and that 
||^llL°°(R d ) < 1- Further, we use the shorthand = (|iV d_2 -UAr) 2 , and assume that 
a>N > 1, for all iV > 1 (this possibly involves the replacement of finitely many Un). We 
then see that for \z\ < Cq 1 and N > 1, 

I[e*~~'">] = K[e«p V) - Hf JV«<V,l>j}] 

(3 ' 16) ™expfciV« E £ {V ,(G N Vr-H)\ 



By Proposition 3.1 and the fact that ^f- N d 2 /a% equals | when n = 2 and is bounded 
by I and converges to zero for N tending to infinity, when n > 3, we find that 



(3.17) IV^] — ► for |*| < c^ 1 

iV-s>oo 



The claim (3.15) then follows by similar arguments as below (3.13). □ 
Remark 3.4. By (1.12) one sees that for any N > 1 



E 



E (i + \y\y (d+1) L Ny , 

2/eLiv 



UN 



< 00. 



Thus on a set of full P-measure, (1 + \y\)~^ d+1 \ and hence all functions V £ S(R. d ), are 
/^-integrable. Redefining C N on a negligible set as (for instance) being equal to zero, we 
see that C N (u) is a tempered distribution for each u£0, and (£ N , V) a random variable 
on (0, A, P) for each V G S{R d ). Hence, we can also view C N as an l S / (M d )-valued random 
variable on (H, A,W). □ 



In the next result about convergence in distribution of C N , the space S'(R. d ) is tacitly 
endowed with the strong topology (see p. 60 of [13], or p. 5, 6 of [17]). We recall the 
notation P G for the law on iS'(M d ) of the massless Gaussian free field (see (1.15)). The 
next result states the convergence of Cn in distribution to the massless Gaussian free field 
in the high intensity regime. 

Corollary 3.5. (under (0.4)) 

(3.18) jC n converges in law to P G as N — > 00. 

Proof. By the Levy continuity theorem on S'(M. d ) (see Theoreme 2, p. 516 of [26], and 
also Theoreme 3.6.5, p. 69 of [12]), it suffices to show that for each V G S(M. d ), the vari- 
ables (C N , V) converge in distribution to ($, V) (under P G ), namely a centered Gaussian 
variable with variance E(V,V), (see (1.15)). 

We thus consider V G 5(M d ) and a continuous [0, 1]- valued compactly supported 
function xl on R d equal to 1 on C 2 (see (3.2) for notation). We write Xl(-) = x{j) f° r 
L > 1 and define Vl = XlV. 

By dominated convergence E(Vl,Vl) — > E(V,V), as L — > 00. By Theorem 3.3 our 
claim (3.18) will thus follow once we show that 

(3.19) lim sup E\(C N , V - V L f] = 0. 

L^oo N>1 
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By differentiating, (3.16) twice in z at the origin and approximating V by Vl>, with 
L' — > oo, we see from Fatou's Lemma that 

(3.20) E[(C N , V - V L ) 2 ] < (V - V L , G N (V - V l )) Ln , for N > 1, L > 1. 

Applying (3.3) to %(•) |V|(- + a) for a G Z d and using that sup{g7v(2/)j V W, |y|<x> > 1} 
is uniformly bounded in N (see (1.7)), we see that p = sup Ar>1 \\Gn \V\ \\l°°(l, n ) < oo, and 
that 

(V-V L ,G N (V-V L )) <cp £ sup \(l- XL (y))V(y)\ — ► 0. 

aez d yea+[o,i) d L ^°° 

The claim (3.19) follows, and this completes the proof of Corollary 3.5. □ 



4 Scaling limits via the isomorphism theorem 

In this short section we revisit Theorem 3.3 under the perspective of the isomorphism 
theorem stated in (0.6). This offers a different route to Theorem 3.3 for sequences (un)n>i 
in a "sufficiently high intensity regime" . The special role of dimension 3 where we recover 
the full range (0.4) will be highlighted. Notation for the Gaussian free field on 7L d have 
been introduced at the end of Section 1. In particular, ((p x ) x ez d stands for the canonical 
field on M zd and P 9 for the canonical law of the Gaussian free field. 

We introduce the random signed measure on M. d 

Given a continuous compactly supported function V on M d and a positive sequence 
(u N ) N >i the identity (0.6) implies that under P 9 ® P 



\ E v(£)(<£-g(0))+ E v(^-)(L x , UN -u N ) 

i^)W-«(«))+VHJEv(i) ft . 



law 



With the notation (0.5) and (4.1), this identity can be rewritten under the form 
E V 9(0)) +V2d^ N^(C N ,V) 



1 \ — > T T ( X \ I •) /r\\\ , /TT1 t>t^4-1 / t r\ law 

2 

(4 ' 3) \ EV^)( V l-g(0))+VM^N^(^,V). 



Importantly, note that C N involves un in its definition, but <& N does not (cf. (4.1)). 

In the next lemma we look at the size of the terms in (4.3) which involve the Gaussian 
free field. To this end, we introduce the sequence 

N 4 when d — 3, 

(4.4) b N = { N 4 \ogN when d = 4, 

N d when d > 5. 
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Lemma 4.1. 

(4.5) ,V) is a centered Gaussian variable with variance (V,GnV) l ^. 

(4.6) (<§> N ,V} converges in distribution to a centered Gaussian variable 
with variance J V(y) G{y — y') V(y') dydy', as N — > oo. 



(47) ^1(E>ti)K-#) 



has a positive limit when N — > oo. 



Proof. We begin with (4.5) and note that (& N ,V) is a centered Gaussian variable with 
variance E* x>&iA V (§ ) 9 ( x , x ') v (jf) (1 '^=*' 6) (V,G n V) Ln , whence the claim. 

(1.8) 

Next, (4.6) is an immediate consequence of (4.5) and (3.4) (with n = 2). We then 
turn to (4.7). We note that from Lemma 5.2.6, p. 201 of [24], one has the identity 

(4.8) E p9 [(<pl - g(0)) (<pl, - g(0))} = 2g 2 (x, x% for x, x' E Z d . 

As a result we find that 

« 4 - 9 » i EP "[{ e/(i^- 9 (o))) 2 ] =£ 

When d = 3, G(y) = c\y\~ l (see (1.3)) is locally square integrable on M 3 , and the above 
quantity equals (see (1.6)) 

(4.10) ^ E V{y)gl{y-y>) V{y>) — ► 18 / V{y) G 2 (y - y') V(y') dydy', 

as can be seen by separately considering the terms \y — y'\ > 7 and \y — y'\ < 7, with 
y,y' G Ltv, in the spirit of what was done below (3.8), letting first iV go to infinity, and 
then 7 go to zero. Note that the integral in (4.10) is equal to (with p t (-, •) the Brownian 
transition density, see above (2.1)): 

2 



dsdtdzdz'i / dyV(y)ps(y,z)pt(y,z')) 



by a similar calculation as in Proposition 4.8, p. 75 of [36]. This quantity is positive when 
V is not identically equal to zero, see the bottom of p. 75 of [36] for a similar argument. 

When d = 4, on = N 4 \ogN, and given L > 1 and 7 > 0, the contribution in the 
right-hand side of (4.9) of the terms with \x — x'\ < L or \x — x'\ > 7^ is O(-^j^), as 
N tends to infinity, by (1.3). Combined with the fact that Eo<|rr|<Ar ^ tends to 

a positive limit, when N tends to infinity (as can be seen by looking at the difference of 
the sum and the integral /e( jv)\B(o 1) T-cf*)' we ^ n< ^ that wnen d = A, 

( 4 - n ) AM1 2 AT E v(^\g\x-^)v(^\ — >c [ V\y)dy. 

Finally, when d > 5, the sum Ea;ez d 9 2 ( x ) converges by (1.3), and hence 

(4.12) A E v(-) g 2 (x-x')v(-) [ V 2 (y)dy. 

We have thus shown (4.7), and this concludes the proof of Lemma 4.1. □ 
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We will now conclude this short section with some remarks, in particular linking 
together (4.3), Lemma 4.1, and Theorem 3.3. 

Remark 4.2. 

1) By similar arguments as in Remark 3.4 we can view Q N as an 5' (Revalued random 
variable. The same proof as in Corollary 3.5 yields that 

(4.13) & N converges in law to P G , as N — > oo. 



2) When u N N d+2 /b N -» oo, we can divide both members of (4.3) by y/2du N N^ +1 . By 
Lemma 4.1 we then see that for each continuous compactly supported function V on M. d , 
both (jC n , V) and (<& N , V) converge in distribution to a centered Gaussian variable with 
variance (V,GV) (once again, note that enters the definition of C N but not that of 
$ Ar ). One thus recovers the convergence asserted in Theorem 3.3 under the assumption 

fiiN N — y oo, when d — 3, 
u n ; — 77 -> oo, when d — 4, 
logiV ' 

un N 2 — > oo, when d > 5. 

This recovers the full range of validity of Theorem 3.3 (i.e. un N d ~ 2 — > oo), when d — 3, 
but only part of the range when d > 4. Observe that when d > 4, in the high intensity 
regime un N d ~ 2 — > oo, (£ iV , V) and ($ Ar , V) have the same distributional limit, regardless 
of whether (4.14) breaks down or not. 

3) One may wonder from the above discussion, whether some cancellations may be ex- 
tracted from (4.3). In the high intensity regime N d ~ 2 — > oo, given M > 1, can one for 
each iV couple two Gaussian free fields ((p x ) xeZ d, {jp x ) x & d i with (L XiUN ) xeZ d, so that 

( ( Px)xei. d is independent of (L x , UN ) x& d, 

- ip 2 x + L XtUN = - (ip x + v / 2mjv) , a.s., for x E NC M (see (3.2) for notation), 



and for any continuous function V supported in Cm, 

\ E V (|) (^ - $1)1 ( ^ iVi +1 ) in probability? 

Such a coupling would yield a simple way to recover the identity of the distributional 
limits of (C N , V) and ($ Ar , V) in the high intensity regime. Of course, when (4.14) holds 
such a coupling can be achieved (by (4.7) the last condition is automatically satisfied). □ 



5 The special case of dimension 3 

We have seen in the last section that when d = 3 the isomorphism theorem (0.6) offers 
a quick route to the study of the asymptotic behaviour of C N in the full range of the 
high intensity regime (0.4). In this section we investigate what happens "at the edge", 
in the constant intensity regime (0.2). The scaling limit of the distributional identity 
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(0.6) (with proper counter terms) will bring as a by-product an isomorphism theorem 
for 3-dimensional Brownian interlacements, see Theorem 5.1 and Corollary 5.3. This 
last corollary has a similar flavor to the distributional identity derived by Le Jan in the 
context of a Poisson gas of Brownian loops at half-integer intensity, see Chapter 10 §2, 
p. 104 of [23]. We recall the notation of the end of Section 1 concerning Gaussian free 
fields. Throughout this section d = 3. 

We denote by H the Gaussian space, which is the L 2 (P G )-closure of {($,/); / G 
S(R d )}. For Y, Z in H the Wick product, see [18], p. 23, 24, :YZ: equals YZ-E pG [YZ], 
and one has the identity, see [18], p. 11, 12 or [24], p. 201: 

E pG [:Y 2 : :Z 2 :] = 2E P [YZ} 2 , for Y,ZeH 

' ' ' = 2(f, Gh) 2 , when Y = ($, /}, Z= <$, h), with /, h e S(R d ). 

When d — 3, the Green function G(y,y') is locally square integrable, and for bounded 
measurable functions V on M 3 vanishing outside a compact set, one can define 

(5.2) J V(y) :<S> 2 y . dy = lim J V(y) :$J >e : dy in L 2 (P G ), 

where <& ye = (*S>,p yi£ }, with p 2/ , £ (-) = p-p(-^), and p(-) a non-negative, smooth, com- 
pactly supported function on R 3 such that f p(z) dz = 1. This fact uses (5.1), the crucial 
local square integrability of G(-), and similar arguments as in the proof of Proposition 
8.5.1, p. 153 of [15], see also [23], p. 101. In addition, the limit object in the left-hand side 
of (5.1) does not depend on the specific choice of p(-). For convenience, we will assume 
that p{ ) is radially symmetric and supported in -8(0, 1). 

For V as above, one also defines (in a simpler fashion) the element of H: 

(5.3) J V(y) %dy = lim J V(y)%, £ dy in L 2 (P G ). 

We now assume that we are in the constant intensity regime (0.2), and for V a contin- 
uous compactly supported function on M 3 , we rewrite (4.2) as the following identity in 
distribution under P 9 ® P: 

(5.4) \ E v(^): V 2 x : + dN 2 (C\V) l = \ E V (£) :(<p x + , 

where :((p x + y/2u N ) 2 : = :ip 2 x : + 2y/2u N tp x + 2u N and :>p 2 x . =<p 2 x - g(0). 

We already know that (C N , V) converges in distribution to (C a , V) (see Theorem 3.2). 
The limit behavior of the other terms in (5.4) is described by 

Theorem 5.1. (d = 3, V continuous with compact support) 
Fora>0 and u N = 3a/N (= da/N d ~ 2 ), as N -> oo, 

( 5 ' 5 ' ) 3^ £), V (jj) : ( LPx + ^ 2u n) 2 '- converges in law to J V(y) :(% + V2a) 2 : dy 

(the last term can be defined as f V(y) :§ 2 : dy + 2y / a J V(y) $ y dy + 2a J V(y) dy, but 
see also (5.10) below). 

Under P G ® P one has the distributional identity 
(5.6) \ J V(y):^ 2 y :dy + (C a ,V) l = 1 - J V(y) : + v^ 2 : dy. 
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Proof. We first observe that (5.5), (5.4) (choosing a = to handle the first term of (5.4)) 
and Theorem 3.2 readily imply (5.6). Thus, we only need to prove (5.5). This is a case 
of lattice field approximation (see chapter 9 §5 and §6 of [15], also Chapter 8 of [33]). 
However, the fact that we consider the d = 3 situation and massless free fields, makes 
the set-up a bit different. Since some care is needed, see for instance [15], p. 185-187, we 
sketch the proof for the reader's convenience. We define 



(5.7) 



<P v ,n = \ Vn v (with d = 3), for y e h N , 



and note that the first term of (5.5) equals 



(5.8) 



^3 E V(y):(<p VtN + y/2Q'' 



We now choose e(N) > tending to not too fast so that 

(5.9) e(N) ->■ and Ne 3 (N) ->■ oo, &s N —■ oo. 

As we explain below, the claim (5.5) will follow once we establish the following three facts: 

(5.10) lim| J V(y):{% + V2^) 2 :dy- J V(y) V^) 2 : */|| L2(pG) : 

(5.11) lim fv(y) :{%, £ + V2^) 2 :dy- -L £ V{y) : {%, £ + ^) 2 : dy 
and for small real z, i.e. \z\ < r , 

¥ (^ G [ exp {^ EV(y):(%, £ + V^) 2 :}] 

N v UV y£h N JJ 

ex p{4 E V^l^ + v 7 ^) 2 :}!) =0 



0, 



Li(P G ) 



= o, 



(5.12) 



- £ 



We first explain how these three facts yield a proof of (5.5). By Theorem 3.50, p. 39 of 
[18], we can replace L l (P G ) by L 2 (P G ) in (5.11). Together with (5.10) we find that 



(5.13) lim -L £ V(y):(%, £ + V2^) 2 :- j V {y) :{% + V2a) 2 : dy 
By Theorem 6.7, p. 82 of [18], we thus see that for some r\ > 0, 



L 2 (P G ) 



0. 



(5.14) 



sup E 

N>1 



pG 



cosh(^ E V(y) + v^) 2 : 



< OO. 



By (5.12), we then find that for some < r 2 < r\ A r , 



(5.15) 



sup £ 

N>1 



Pa 



cosh(^ E ^(vVv + v^ 2 :)" 
vv yew /J 



< oo. 



Thus, by a similar argument as below (3.12), we see that the laws under P 9 of the random 
variables -^3 EyeLjv V(v) '•( L Py,N + v^a) 2 : are tight, and, when along a subsequence 7V fe 
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they converge in distribution to a random variable U, then, for \z\ < r 2 , e zU is integrable 
and one has 



zU 



(5.16) 



1 = lim E pa ' 


N k 




(5 i 2) hm E P ° 






= E pG 


exp 



cxp 



z 

7V3 



E V(y):(ip y , N + V2^) 2 :} 

y&L N J 



using in the last step both the uniform integrability of the variables under the expectation 
in the second line when \z\ < r 2 by (5.14), and (5.13) (see Theorem 5.4, p. 32 of [3]). 
By analyticity, the equality between the first and the last term of (5.16) extends to the 
strip {z G C; \Rez\ < r 2 }. Hence, the random variable U has same distribution as 
JV(y) :(& y + \f2a) 2 \ dy, regardless of the extracted subsequence N k , and (5.5) follows. 

There remains to establish (5.10) - (5.12). The claim (5.10) is an immediate conse- 



quence of (5.2), (5.3) and the identity :($ y , e + V2a) 2 : = :® 2 y£ : + 2y/2a $ y , £ + 2a. 
We then turn to (5.11) and simply explain why 



(5.17) 



lim 

N 



( V{y):<S>l e :dy-± £ V(y):& w 

J y£h N 



The case of the linear term (where 2v / 2a <& y £ replaces :^ £ is simpler to handle, and 
the constant terms (after developing the square in (5.11)) cancel each other. The quantity 
under the limit in (5.17) is bounded above by I x + 1 2 , where 



(5.18) 



h = E 1 



(J h N {y)^l £ :dy) 



, with 

G z + 



Ml/)= E {V(y)-V(z))l{yez + ±-iO,ir}, 
zeh N k JV J 

e- ~ -Ke'- d y\]- 



h= E \V(z)\E 



^+^[0,1) 3 



(5.19) 
By (5.1) we see that 

(5.20) I\ < 2 J \h N {y)\ \h N {y)\ G 2 £ (y - y') dydy', where 

(5.21) G £ (y) = E pG [^ £ $ 0)£ ] = f Pe(z - y) G(z - z') p £ (z') dzdz', for y G M 3 . 

Note that p £ (-) is spherically symmetric with support in B(0, e) and G harmonic on 
M 3 \{0}, so that 



(5.22) 



G e {y) = G(y), when \y\ > 2e 



< 



3 , when \y\ < 2e, 



where in the last step we used that f\ z _ z ,\ <r G(z — z') dz = cr 2 , for r > 0, z' G R 3 . 
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We introduce a box A = [— M, M) 3 , with M > 1 integer, such that 

(5.23) A contains all points within | • loo-distance 1 from the support of V. 

Note that vanishes outside A (cf. (5.18)), and converges uniformly to as A" — > oo. 
Moreover, by (5.22) and the local square integrability of G 2 , we see that j^ x ^G 2 {y — 
y') dydy — > j AxA G 2 (y — y') dydy', as A" — > oo. Hence, see (5.20), 

(5.24) h -> 0, as N -> oo. 

We then control J 2 - We first note that A^ 3 |A| bounds the number of boxes z + -^[0, l) 3 
intersecting the support of V (with |A| the volume of A, and z G h N ). Hence, 

(5.25) J 2 < llVHco |A| sup E pG [\:§l £ : - :^ £ :|]. 

j/e^[o,i) 3 

Moreover, for y G M 3 we have 

E pG \^l £ : - :^-f] = E pG [.^l £ :] + E pG [:$l:] -2E pG [:^ £ : 

(5.26) ^ 4(G £ 2 (0) - G 2 (y)) { T ^(G e (0) - G £ (y)) 

< C M^ (since IIIV^HU^s)^^). 

We thus see that 

(5.27) I 2 < c(p) \\V\UA\ (iV5 3 )"l — » (by (5.9)). 



N 



This complete the proof of (5.11). We now turn to (5.12). 

By (5.22) and (1.2), (1.3), (1.6) we have for y,y' Eh N ,N> 1, 

G e {y-y') <c(p)(\y-y'\ V -M 'and 

(5.28) V 

9N(y,y')<c(\y-y'\V-) . 

We can thus make sure that for small real z, 

(5.29) sup M £ G £ {y-y') \V{y')\ < \, and 

(5.30) sup M E Wy-y')IW)l<^ 

(incidentally, by (3.6) and the continuity of G \V\, we also have \z\ \\G \V\ ||loc( M 3) < | 
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If An = AflLjv, the matrix (G £ (y — y')) y , y 'eA N is positive definite. Indeed, by (5.21) it 
suffices to note that p £ (- — y), y G A^r, are linearly independent (one easily sees by taking 
Fourier transforms and using an analytic continuation argument that any linear relation 
between these functions is trivial). 

By Lemma 5.2.1 on p. 198 of [24], or Proposition 2.14 on p. 47 of [41], for z as above, 



(5.31) £ pG [exp{-^3 £ V(y) :(%, £ + 



^{V^I-zGeV)- 1 !)^ 



(det 2 ((J - zG £ V)\ KnxKn ))* 



i •> 



(5.32) ^ P9 [ex P {^3 + 



(det 2 ((I - zG N V) lANxAN )y- 



where det 2 ((/ — zG £ V)\a n xA n ) denotes the regularized determinant of the finite ma- 
trix K, which is the restriction to A^ x A^r of {8 y , y > — zG £ (y — y')^w^) yy 'eL N -> * s 
(det K)e~ Tr ( K ~^ (where Tr stands for the trace), see [34], p. 75. In addition, with some 
abuse of notation, in the term in the exponential in the right-hand side of (5.31), G £ 
stands for the linear operator on functions on h N defined by a similar formulas as (1.8), 
with <7jv(2/j y') replaced by G £ (y — y'), so that G £ V operates on bounded functions on h N 
and I — zG £ V is invertible by (5.29). Also, compared to the above references, we used 
the identity following from the Neumann expansion of (I — zG £ V)~ l : 

(V, (I - zG £ V)-H) Ln = (V, 1) Ln + z(V, (I - zG £ V)-'G £ V l)^. 

In a similar manner, with analogous notation, for z as above, we have 

V(y):{% !N + V2^Y:\ =7— 

By Proposition 3.1 and (5.30) we have for z as above 

(5.33) lim (V, (/ - zG N V)-H) hN = (V, (I - zGVyH). 

By an analogous argument as for the proof of Proposition 3.1 and (5.29) we have 

(5.34) hm (V, (/ - zG £ V)- l l\ N = (V, (/ - zGV)' 1 !). 
As we now explain, the claim (5.12) will now follow once we prove that 

(5.35) sup £ (G £ (y, y'f + g N (y, y'f) 4 < oo, 

N y,y> & A N ^ 

(5.36) hm £ {Gc(y,y')-g N (y,y')) 2 ^ = 0. 

N y,y'eA N iV 

Indeed from Theorem 9.2, on p. 75 of [35], one knows that for any nx n matrices A, B 
(and || • || 2 denotes the Hilbert-Schmidt norm) 

(5.37) |det 2 (/ + A)-det 2 (/ + S)| < \\A - B\\ 2 exp{c(||A|| 2 + \\B\\ 2 + l) 2 } 
(the constant c does not depend on n). 

Assuming (5.35), (5.36) one thus finds that for small z 

(5.38) lim det 2 ((J - zG £ V)\a nxAn ) - det 2 ((J - zG n V)\ AnxAn ) = 0, and 

(5.39) the regularized determinants in (5.38) each belong to [|, |] for N > 1 
(we use (5.37) with B = and det 2 (7) = 1 for this last fact). 
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Combined with (5.33), (5.34) the claim (5.12) then follows. 

There remains to prove (5.35) and (5.36). The bound (5.35) follows from (5.28) in a 
straightforward fashion, and we concentrate on (5.36). We pick 7 > 0. The sum in (5.36) 
equals J\ + J 2 , where J\ collect the terms where \y — y'\ > 7 and I2 the terms where 
\y — y'\ < 7- When N is large, G s (y — y') = G(y — y') for all \y — y'\ > 7 (see (5.22)), 
G(y — y') — QNiy — y') converges uniformly to zero on this set by (1.7). We thus find that 
lim^ Ji = 0. As for the second term J 2 , by (5.28) we have 

(5.40) hW 2 <IS%^ (n 2 + £ *)<J(A,p) r 

Letting 7 tend to 0, we obtain (5.36). This completes the proof of (5.12) and hence of 
Theorem 5.1. □ 

Remark 5.2. 

1) Note that for V continuous with compact support on R 3 one has 



E 



2n 



V(y):& y :dy) \=2j V(y)G 2 (y - y')V{y') dy dy' < 

l^ / l|L 1 (R d )(||^ / ||L°°(R d ) + II^IIlMK 3 ))' 

and a similar bound holds for E[( J (V(y) $ y dy) 2 ] = J V(y) G(y — y') V(y') dy dy'. 

We can thus extend the linear map sending V to J V(y) :(<f> y + ^2a) : dy e L 2 (P G ) 
by continuity to 5(R 3 ) for any a > 0. By the regularization theorem, see Theorem 2.3.2, 
p. 24 of [17], one can find a version denoted by :($ + V2a) 2 :, which is an <S'(IR 3 )-valued 
random variable defined on the canonical space where P G is defined such that for any 
V e S{R 3 ), 



> G -a.s., <:($ + v / 2^) 2 :,V} = J V(y) :(% + V^f: 



dy. 



2) The intensity measure of C a equals ady (see (2.39)), and by similar considerations as 
in Remark 3.4 we can also view C a (after a possible redefinition on a set of measure 0) as 
an S' (R 3 ) -valued random variable defined on (f2, A, P). □ 

As a consequence of Theorem 5.1 we obtain an isomorphism theorem for the occupation 
measure of Brownian interlacements on R 3 , which arises as the scaling limit (with suitable 
counter terms) of (0.6): 

Corollary 5.3. (d = 3,a>0) 

Under P G ® P one has the identity in law on iS'(R 3 ): 



(5.41) \^ 2 : + C a l =\:^ + ^) 2 - 



Proof. By Remark 5.2 we can obtain the identity (5.6) for any V G 5(R 3 ). Hence under 
P G <g) P, the random variables (| :$ 2 : + C a , V) and (| :($ + \f2a) 2 \, V) have the same 
distribution for each V G <S(R 3 ), and (5.41) easily follows. □ 

Remark 5.4. In the context of Brownian loops on R d , 1 < d < 3 (suitably killed when 
d = 1,2), Le Jan, see p. 104 of [23], showed that | :$ 2 : (or its corresponding object when 
d = 1,2) has the same law as the renormalized occupation field of a Poisson cloud of 
Brownian loops with intensity parameter ol = \. 

Observe that here, in (5.41), no renormalization on C a is needed (and a varies over 
R + ). Only terms involving the free field need renormalization. □ 
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